Gapless superconductivity and the Fermi arc in the cuprates 
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We argue that the Fermi arc observed in angle resolved photoemission measurements in under- 
doped cuprates can be understood as a consequence of inelastic scattering in a d-wave supercon- 
ductor. We analyze this phenomenon in the context of strong coupling Eliashberg theory, deriving 
a 'single lifetime' model for describing the temperature evolution of the spectral gap as measured 
by single particle probes such as photoemission and tunneling. 

PACS numbers: 74.20.-z, 74.25.Jb, 74.72.-h 



Underdoped cuprates are characterized by a gap which 
exists in the excitation spectrum even above the super- 
conducting transition temperature T c and is detectable 
up to a temperature T* , the "pseudogap temperature". 
The origin of this pseudogap has been a source of much 
debate. One of the most intriguing experimental charac- 
teristics of this pseudogap regime is a partial destruction 
of the normal state Fermi surface: states near the edges 
of the two-dimensional Brillouin zone ('antinodes') are 
gapped below T* , while states along an arc near the zone 
diagonals ('nodes') remain gapless. This phenomenon 
has been extensively measured in angle resolved photoe- 
mission (ARPES) experiments 1 ^! 3 -. The convention used 
to interpret ARPES data is that a point \tp on the Fermi 
surface is considered to be gapped if the spectral function 
A(k.F,uj) cx |ImG(kf,w)| has a maximum at a non-zero 
frequency and is gapless if the maximum of A(kp,o;) is 
located at zero frequency^. In underdoped cuprates in 
the temperature range between T c and T*, the gapless 
regions appear to form arcs rather than closed contours. 
The most recent measurements show 3 that the length of 
the arc is temperature dependent: it appears to be linear 
in T and seems to vanish if the data obtained in the tem- 
perature range T c < T < T* are extrapolated to T = 0. 
Moreover, as T is decreased below T c , the arc collapses 
rapidly to a point node^. Understanding these Fermi arcs 
has become a key question in the physics of the cuprates. 




Theoretical scenarios for the arc can be separated into 
two groups. One scenario is that the arc emerges due 
to long range (or quasi-long range) ordering that is not 
related to the pairing. The proposals include orbital cur- 
rents^, SDW&l, CDW£, and dDW-2 ordering. A generic 
consequence of these theories is the presence of two gaps 
- a smaller superconducting gap in the nodal regions, and 
a larger gap of different origin in the antinodal regions. 
Above T c , the superconducting gap vanishes leaving only 
the other order, which leads to Fermi surface pockets near 
the zone diagonals 10 . The ARPES intensity on one side 
of these pockets is expected to be weaker than on the 
other side (as observed in experiments on materials with 
density wave ordering) so that the pockets look like arcs 
in ARPES measurements. The notion of two gaps is also 
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FIG. 1: (a) Im G around the Fermi surface from Eq. ^ with 
$ = $o cos(2<^)), where $o=40 meV and 7 = 30meV. The 
top curve is at the node (<j> = 45°), the bottom curve at the 
antinode (<f> = 0°), with cj> increments of 5°. (b) Spectral gap 
(half the peak to peak separation) versus cj> for 7 ranging from 
(dashed curve) to 60 meV (bottom curve) with 7 increments 
of 10 meV. 



implicit in scenarios relating the arc to Mott physic s 11 ' 12 . 

An alternative 'one-gap' scenario interprets the pseu- 
dogap as a remnant of the superconducting gap, with 
thermal fluctuations destroying long-range order at T > 
T c , but allowing gap- like features to remain in the spec- 
tra for a range of higher temperature s 13 ' 14 ' 15 . The effect 
of phase disordering on the spectral function has been 
extensively discussed in Refs. I2lfl6l A phenomenological 
Green's function has the form 



G(k 



1 



uj + 271 



«72 



W-M71 



i 7 i)(w -M72) - $ 2 (0) 

Here $(</>) is the d— wave anomalous self energy, and <p 
is an angle denoting the position along the Fermi sur- 
face. The scattering 71 comes from pairing (or vortex) 
fluctuations, and 72 is the fermionic scattering rate. 

A simplified model with 71 = 72 = 7 is sufficient to 
understand how the Fermi arc, as experimentally defined, 
emerges due to scattering. At small 7, ImG(kp,w) has 
sharp peaks at u> — ±$(</>). However, as 7 increases, 
the peaks broaden and slowly shift to smaller absolute 
frequencies (Fig. la). For 7 < -\/3| < I ) (</ , )|, the maxima 
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are located at w = ±[2<J>a/7 2 + $ 2 - -f 2 - $ 2 ] 1 / 2 . For 
7 > \/3|<j>(0)|, a single maximum occurs at lo = 0. The 
boundary value |$(</>)| = j/VS sets the length of the 
Fermi arc (Fig. lb). Because $(0) vanishes at the node, 
the arc has a non-zero length for arbitrarily small 7. 

Our goal is to provide a microscopic justification for 
Eq. |T]) (with 71 = 72 = 7) and to show that this is 
equivalent to gapless superconductivity. The latter nat- 
urally emerges at non-zero T in an Eliashbcrg theory, 
independent of the nature of the electron-boson interac- 
tion and the symmetry of the pairing gap. Obtaining 
a two-lifetime model with 71 ^ 72 would require addi- 
tional physics£ii£. We show that even in the absence 
of impurities, 7 is always non-zero for T > 0. Both 7 
and the length of the arc scale linearly with T at high 
temperatures, when the strongest interaction is between 
fermions and classical (thermal) bosons. The linear T 
behavior of 7 sets in well before the temperature exceeds 
the "effective Debye frequency" of the bosons for the 
same reason that linear T resistivity becomes apparent 
for T > ft Debye /3i£. The linear T dependence of 7 in the 
pseudogap phase could also arise from the same physics 
that gives rise to the marginal Fermi liquid behavior of 
the normal stated. 

This scenario is consistent with the ARPES experi- 
ments. The key evidence is the vanishing of the extrap- 
olated arc length at T = 0, and the linear T dependence 
of the arc length for T > T c . Consistent with this is the 
"closing" of the gap in the nodal region, i.e. the reduction 
of the spectral peak position to zero for small $ relative 
to 7, and the "filling in" of the gap in the antinodal re- 
gion, i.e. the spectral peaks broadening without moving 
much for large <£> relative to 7. This behavior is evident 
from Eq. fl} (see Fig. la). For a two-gap scenario, the 
arc length above T c should have only a weak dependence 
on T, and would extrapolate to a non-zero value in the 
T=0 limit, in contradiction to the ARPES datai -. 

The rapid shrinking of the arc length below T c found in 
Ref. is a more subtle issue. Our analysis shows that a 
Fermi arc with T-dependent length is a generic property 
of a d-wave superconductor at non-zero T, independent 
of issues associated with superconducting coherence. In 
other words, long-range superconducting order does not 
have to be destroyed by fluctuations for the arc to be 
present. This generally implies that the arcs must survive 
below T c . We find, however, that the length of the arc 
is very small at the lowest T (exponentially small if the 
spectrum of the scattering bosons is completely gapped 
well below T c ), but becomes both sizable and linear in 
T at temperatures above the fluctuation-driven T c . This 
would be consistent with the drop of the scattering rate 
below T c that has been inferred from a variety of tech- 
niques, including ARPES and various conductivity mea- 
surements^. Another possibility, which we don't explore 
here, is that the rapid variation of the arc length near T c 
is due to a strong temperature variation of the damping 
originating from the interaction between electrons and 
pairing and/or vortex fluctuations^ ' 16 ' 20 ' 21 . Phenomeno- 



logically this can be described as the vanishing of 71 in 
the two-lifetime model of Eq. fl]). Once 71 = at T c , 
A(]s.f,u>) vanishes at lo — 0, and the spectral peak is lo- 
cated at a non-zero frequency for all kp points except 
the d-wave node^. 

Gapless superconductivity was first discussed in the 
context of a BCS s-wave superconductor with magnetic 
impurities^. Although not widely appreciated, even in 
a clean s-wave system, the Eliashberg equations yield 
gapless superconductivit y 24 ' 25 ' 26 ' 27 . The reason is that 
at non-zero T, the density of the bosons is non-zero, 
and a fermion can undergo scattering by on-shell ther- 
mal bosons 26 . This leads to a non-zero damping term in 
the fermionic self-energy £(<*>), just as in a superconduc- 
tor with magnetic impurities. 

To understand gapless superconductivity for the clean 
case and to discuss the Fermi arcs, we consider a d-wave 
Eliashberg theory with an effective (dimensionless) in- 
teraction x(k — k', fl) between the electrons at the Fermi 
surface (|k| = |k'| = kp). We assume x(k — k',fi) can be 
decoupled into s-wave and <i-wave harmonics, i.e., 

X(k - k', fl) — X [xd{£l)d k d k > + Xs{tyskSk>] (2) 

where A is a dimensionless coupling, Sfc = 1, and d k = 
i (cos fez — cos k y ). We choose a normalization such that 
Xs(0) = 1. We assume that \d is attractive, and Xs 
is repulsive. We also assume for simplicity that there 
is no feedback effect from the pairing on x(k — k',f2) - 
this last simplification may be partly justified by the fact 
that, experimentally, feedback effects are relatively weak 
in the pseudogap phase. 

Eliashberg theory provides a mechanism for the open- 
ing of a gap in the excitation spectrum, without giving 
information on phase coherence. Therefore, in what fol- 
lows, we shall use it to describe the pseudogap behavior 
for T < T*. In a superconductor, the self energy has 
a normal part expressed in terms of a renormalization 
factor Zk(oj) (which is k independent if Sk — 1), and an 
anomalous (pairing) part expressed in terms of the am- 
plitude $fe(w) = dk$(ui) (the equality follows from our 
choice of interaction, Eq. [2]). The superconducting gap is 
A fc (w) = <&fe(w)/Z fe (w) = d k A(io). The quantities A(lo) 
and Z(lo) are given by the solution of two coupled in- 
tegral equations^, of which the equation for A involves 
only A while the other expresses Z in terms of A. As 
long as the interaction is finite at zero frequency, $(cj m ) 
tends to a non-zero value at the smallest u> m , and for the 
purposes of our small lo analysis can be safely replaced 
by a frequency- independent value <E>. The low- frequency 
behavior of Z(L0 m ) is more involved, and the conversion 
of Z to real frequencies at finite T is nontrivial 27 . 

Note that Eq. {2} is an approximation - there is no 
guarantee that for a realistic x(k — k', f2) the pairing am- 
plitude $ will have the simple d— wave form $(/c) = d k &- 
In particular, near the onset of pairing, the form of the 
gap will be such as to minimize the pairbreaking effects 
specific to systems with a non-s-wave gap symmetry^ 5 -. 
This by itself may give rise to a suppressed gap around 
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the node, even without the contribution from Z^{uj). We, 
however, will not explore this possibility here and simply 
assume $fc oc d^. 

It is more convenient to solve the gap equation directly 
along the real frequency axis and to cast the Eliash- 
berg equations in terms of D(u>) — A(u))/u>. For an 
electron-phonon superconductor, this procedure has been 
exploited in Ref. 27 to obtain an equation for A(w) along 



the real axis. The extension to a d-wave superconductor, 
and to arbitrary Xs(^) and Xd(ty, is straightforward 28 
and yields 

D(lu)B(uj) = A(u) + C(lu) (3) 

where 
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All integrals should be understood as principal parts. 
The terms A(u) and B(u>) are non-singular at small fre- 
quencies and can be safely approximated by A(u>) — 
A(0) = const, and B(ui) = B\ oj (B\ = 1 + A for vanish- 
ing -D(O)). The non-trivial physics is due to the presence 
of the C{ui) term in Eq. Q. This term originates from 
the imaginary part of the interaction %(k — k', O) along 
the real frequency axis and describes scattering by on- 
shell thermal bosons. Neglecting C(cu) would lead to a 
conventional BCS-type result A(oj) = A(uj)uj/ B(uj) — 
A(0)/Bi = const. However, the second term in the inte- 
gral for C(ui) is proportional to D{lo) — A(w)/w, and 
diverges at zero frequency if A(0) is nonzero. Keep- 
ing only this term and approximating the remainder of 
the integral by its oj — > limit, we find from Eq. Q 
C(u) = -iD{uj)Q{T), where 



Q(T) = A r 

J — t 



dCl Imxs(fi) 
sinh p y/l -D 2 (Cl) 

Upon substituting into Eq. §5§ we get 



AH 



Bi 



I'y 



7 = 



Q(T) 
Bi 



(5) 



(6) 



This is a gapless superconductor - the gap A is imag- 
inary and linear in frequency at small frequencies, but 
recovers to a real value, equal to the pairing amplitude 
$, at frequencies larger than 7. We then obtain at small 
frequencies Z(u) = Zo(l + ry/u/), where Z = B\<§>/Aq is 
a mass-renormalization factor. For the Green's function 
we then get 
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This is precisely the d-wave single-lifetime model of 
Eq. (P) (7! =72). Note that the gap amplitude <fr/Z 
and the scattering rate 7 are determined by different in- 
teractions. In particular, 7^0/$ = jBi/Aq is strongly 



enhanced if the repulsive s-wave component of the inter- 
action is much larger than the d-wave component. 

As stated earlier, the model of Eq. gives rise to arcs 
of the Fermi surface simply because the condition for the 
arc, 7 > \/3($/Z )\dk\, is always satisfied sufficiently 
close to the nodal direction, as long as 7 is non-zero, as 
it is for any non-zero T. 

To understand the temperature dependence of the arc 
length, we note that below the onset of pairing, $ and 
Zq are weakly T-dependent, so that the T dependence 
of the arc length comes primarily from j(T). This T 
dependence is given by Eq. ([5]), and the functional form 
of j(T) depends on the form of X«(0). As an example, 
consider the situation where well below T*, Imx s (fi) can 
be approximated by an Einstein mode at some frequency 
f2o, i.e. Imx s (n) — ttCIq S(n 2 — fig). This can be a phonon 
or collective spin excitations - the latter form a gapless 
continuum in the normal state, but become a mode (plus 
a gapped continuum) once $ is nonzero. 

Substituting ImXs(^) into (UJ) and ([5]), we obtain that 
the boundary condition for the arc reduces to 
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where 77 is the normalized density of states at lu — 
Qq. For a circular Fermi surface, this is r\ = 

(2/iT)Re J* /2 d4>/^l -D 2 ((t>,Q ) fa > 1 if the largest 
contribution to r\ comes from antinodal fermions). The 
quantity <&/Zq is the measured gap at frequencies well 
above 7, and it is well approximated by the maximum 
value of the ARPES gap at the antinode. We see from 
Eq. ((8|) that the arc length is exponentially small at small 
T <C fio, but becomes linear in T at high temperatures. 
This by itself is not a surprising result as at large T, the 
interaction is dominated by the thermal (zero-Matsubara 
frequency) term in which case the damping scales as T 
(the zero-Matsubara frequency contribution corresponds 
to approximating sinh(f2/T) in Eq. ^ by Q/T). A more 
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FIG. 2: (a) Fit of j(T) as given by the left hand side of 
Eq. (|8]) to values extracted from experimental ARPES data 29 . 
A constant term of 2.3 meV has been added to account for 
residual/instrumental broadening at T—0. The fitted value 
of fio is 28 meV. (b) Arc length versus T from the fit in 
(a) assuming ($//7o)dfc = $o cos(20) with $o = 40 meV and 
T* = 3007V Solid points are the data from Ref. S The 
dashed curve assumes instead that 7 oc T for all T. 

interesting result, well known from the electron- phonon 
literature^, is that, numerically, the r.h.s. of Eq. © is 
essentially linear in T already from T ~ 0.3£7o- 

Experimentally, the arc length is linear in T for T c < 
T < T*. This is roughly consistent with Eq. |5]| if we 
associate T c with the mode energy (e.g. for the neutron 
resonance mode 30 , T c ~ 0.2f2rj, in which case a linear 
in T behavior sets in at T > 1.5T C ). Moreover, this 
behavior will also exhibit T/T* scaling as in Ref. y if 
&/Zo scales with T* , as would be expected if T* were 
associated with pair formation. To look into this fur- 
ther, we show in Fig. 2a a plot of j(T) as given by the 



l.h.s. of Eq. {SJ versus the experimental j(T) values ex- 
tracted from ARPES^ on a slightly underdoped sample 
of Bi 2 Sr 2 CaCu 2 8 (T C =90K). A good fit is obtained with 
a reasonable value for fig (28 meV) and a prefactor of 
one (corresponding to 77A / (1 + A) = 1). The resulting 
arc length is shown in Fig. 2b. The agreement is quite 
reasonable given the number of approximations made in 
deriving Eq. {SJ, and the fact that we did not take into 
account any temperature or doping dependence of &/Zq 
and £lo when constructing these plots. We note that a 
dependence of 7 that is strictly linear in T is more consis- 
tent with the data, though a better fit can be obtained if 
dk is flattened near the nodes as in Ref.0. In addition, all 
the data shown in Fig. 2b were restricted to temperatures 
above T c (unlike in Fig. 2a). 

In conclusion, we find that the temperature dependent 
Fermi arcs observed by photoemission in the pseudogap 
phase can be explained by lifetime broadening of a d- 
wave paired state. We have investigated this in detail by 
examining the effects of inelastic scattering within the 
context of strong coupling Eliashberg theory, and find 
that the results are in accord with the data. This suggests 
that T* can be thought of as an energy scale associated 
with d-wave pair formation. 
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